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We have performed a new evaluation of the hadronic contribution to = [g — 2)/2 of the 
muon with explicit correlations of systematic errors among the experimental data on a{e'^e~ 
hadrons). Our result for the lowest order hadronic vacuum polarization contribution is ajj"'' — 
702.6(7.8)(14.0) x 10"^° where the the first error is statistical and the second is systematic. The 
total systematic error contributions from below and above = 1.4 GeV are (13.1) x 10"^" and 
(5.1) X 10"^'^ respectively, and are hence dominated by the low energy region. Therefore, new 
measurements on a{e'^e~ — > hadrons) below 1.4 GeV can significantly reduce the total error on 
a'lf^. In particular, the effect on the total errors of new hypothetical data with 3 % statistical and 
0.5 - 1.0 % systematic errors is presented. 
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A new measurement of the anomalous magnetic mo- 
ment of the muon, = {g—2)/2, to an absolute accuracy 
of aeSp - ± 4.0 X 10-1° is proposed by the AGS E821 
Cohaboration at Brookhaven National Laboratory [0 . 
The theoretical value of the muon g—2 value consists of at 
least the three standard model contributions: quantum 
electrodynamics (QED), electroweak (EW) and hadronic. 
The latter arise from hadronic vacuum polarization ef- 
fects caused by effective photon couplings to hadrons via 
charged quarks and consequent quantum chromodynam- 
ics (QCD) interactions with gluons. 

Any residual difference between the sum of the stan- 
dard model contributions and the new experimental value 
^exp ^^Yl be indicative of new physics: 
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The new experimental value can only be sensitive to 
electroweak and possibly supergravity j|] Q and muon 
substructure effects Q provided the errors on the stan- 
dard model contributions are known better than the ex- 
perimental accuracy. The QED and EW contributions 
have been calculated from theory and are known an or- 
der of magnitude better than the expected experimen- 
tal accuracy {a^^^ = 11658470.6 ± 0.2 x IQ-^^ [§ and 
a^^ = 15.1 ± 0.4 X IQ-i" [§). At the same time, the 
lowest order hadronic contribution cannot be calculated 
accurately enough by QCD and hence a phenomenolog- 
ical procedure must be used for its calculation. While 
the results of the published ajj'"' calculations have errors 
which vary from just above E821 error to greater than 
the electroweak contribution (see Tab. ||) and even differ 
in principle of approach, it is of interest for the interpre- 
tation of a new experimental measurement of the muon 



g — 2 value to investigate in more detail the precise pro- 
cedure used for calculation of the hadronic contribution 
total error. 

Fortunately the hadronic contribution to the muon 
g — 2 value can be related to a dispersion integral over 
the experimental cross section a{e^e~ — > hadrons) (see 
Sect. 0). The error on the hadronic contribution to is 
hence determined by the error on the experimental cross 
section and is by all acounts larger than the expected 
E821 accuracy. Therefore, the hadronic contribution er- 
ror will largely determine the sensitivity level of E821 to 
new physics: 
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This is the motivation for several recent evaluations of 
the lowest order hadronic contribution to which are 
presented in Tab. ||. These calculations may roughly be 
placed into two categories of approach whether they are 
based primarily upon (aggresive) model dependent tech- 
niques which yield relatively smaller errors j|] j|] ||l^ or 
on (conservative) model independent techniques (trape- 
zoidal integration) which yield relatively larger errors |ll| ] 
. A typical method for evaluating systematic error is 
comparison of the model-dependent calculation with one 
based on trapezoidal integration. 

While the model dependent calculations are based on 
theoretical innovations to represent various components 
of the cross section a{e^e~ — > hadrons) they lack the 
merit of a clear prescription for correlating the system- 
atic errors of the experimental input. These errors are 
in some cases as large as 20 % in overall normalization 
of the cross section. Therefore, the criterion used to 
select a set of fit parameters, although indicative of the 
best overall fit to the data, may perhaps not indicate the 
true error on the integral over experimental data points 
which is the hadronic contribution to muon g — 2. In 
particular, a shift of the whole curve up or down may fall 
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within the experimental data point errors and yet yield 
a variation in a^'"' larger than the stated errors. 

The model independent techniques for calculating a'^f"^ 
are based upon trapezoidal integration over the exper- 
imental data points. As will be described below, this 
approach allows for a definite procedure for determining 
both the central value for a!*'"^ and its errors both of 
which make explicit use of the statistical and systematic 
errors for each experimental data point input. In par- 
ticular, for each experiment an error weighted fraction 
will be defined and used in averaging over different ex- 
periments which make measurements of the same cross 
sections over common energy regions. 

The most recent model independent analysis [jl^ was 
the first to account for systematic error correlations 



among the experimental data used in the 



^had 



calcu- 



lation. The systematic errors were correlated at each 
energy point separately in the process of determining 
an error-weighted average cross section ratio R{s) (see 
Sect. H before performing the g — 2 integral. Special 
care was taken to calculate the p-meson contribution sep- 
arately from the rest of the hadronic production cross 
sections and there are many significant contributions to 
the discussion of ajj°'^ calculations in this important com- 
munication. 

The main new aspect of the present ajj'"' calculation 
is that the systematic errors are correlated by a different 
error combination formula a/fer performing the g~2 inte- 
gral. Furthermore, correlations are accounted for across 
energy regions and between experimental measurements 
of the p-meson and of the additional hadronic cross sec- 
tions which contribute to cr(e^e~ — > hadrons). (It ap- 
pears that by treating the p-meson separately the previ- 
ous calculation did not correlate some errors completely.) 

In short, the g — 2 integral has been performed over 
each hadronic cross section separately without combining 
them into their equivalent R{s) value for the present cal- 
culation. The focus is therefore upon the error-weighted 
average of the ajj'^'^ integral value itself (and its errors). 



These innovations have yielded a result with total error 
remarkably equal to that obtained previously |l2[ . How- 
ever, there is a difference obtained as to how much of 
the error derives from the energy region above and below 
1.4 GeV which calls attention to the precise procedure of 
error combination employed in the calculation. 



II. BACKGROUND ON Aff'^'^ CALCULATIONS 



The hadronic contribution to consists of a dominant 
lowest order term, shown in Fig. |l|, several higher order 
terms in Fig. |^ (the number above each diagram indi- 
cates how many contributing diagrams are in its class), 
and finally a group of hadronic light-by-light scattering 
terms. Detailed calculations are given elsewhere of the 
Figure |^ and hadronic light-by-light contributions. Since 
their errors are well below the expected E821 error of 
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-5.2 ± 1.8 X 10 
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-9.0±0.5xl0~i" [JJ and 
I) they are not discussed 



here. There is another hadronic light-by-light calculation 
with a larger error a'^^^ht ^ _9.2±3.2 x lO^^" @ but it 
does not depend on the cross section cr(e+e^ — > hadrons) 
and therefore is not of concern to the present calculation. 
In view of the difference between the two latest hadronic 
light- by- light contributions, and the consequent ambigu- 
ity over defining the total hadronic contribution, this pa- 
per is concerned only with the lowest order hadronic con- 
tribution to the muon g — 2 value shown in Fig. |^. 

A. Formalism of the Hadronic Contribution aj^"'* 

The largest contribution to aj^°'*, shown in Fig. |l|, can 
be related to the total Born cross section (lowest order 
in QED) for hadron production in electron-positron an- 
nihilations, cr^^^ = a^{e~^e^ — s- hadrons), by means of 
dispersion theory and the optical theorem . Defining 
^ = s/to^ and /3 = yl — 4/C the result is 
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where the kernel function K{s) in general arises from a massive photon propagator in the g — 2 Schwinger calculation: 
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The cross section ratio R{s) = (7°^J(J°^ with cr°^ = 
47ra^/3s and the kernel K2{s) is used in the R{s) formu- 
lation of the ajj"'* dispersion integral: 

KM ^ ^A-(.). 

The ahad formulation is useful for low energy data 
which are usually published as individual exclusive 
hadronic cross sections; the R{s) form of the ajj'"' dis- 
persion integral is useful for higher energies where exper- 
imental data are usually published as the inclusive ratio 
R. 

From the ajj'^'* dispersion integrals it is apparent that 
the error on the hadronic cross sections determines the 
error on ajj:'"*. Therefore, the AGS E821 experimental 
error goal just larger than 3.5 x 10~^° on the hadronic 
contribution of roughly 700 x 10~^° requires a 0.5 % ac- 
curacy on the hadronic contribution calculation. 



III. THE WFSA EVALUATION PROCEDURE 

A. Correlation Postulates 

As is implied by the g — 2 dispersion integral (Eq. |^) 
the ajj"'^ calculation procedure consists of some combi- 
nation of the following steps: 1) integration over energy, 
2) weighted average over detectors and 3) sum over ex- 
clusive hadronic modes which contribute to (T°(e+e~ 
hadrons). The sequence of these three sums may be 
interchanged as for example performing the energy in- 
tegration last ||l2|. Most calculations, however, employ 
the sequence (1,2,3) — the most natural one — where in 
the last step, the quadrature sum of errors over modes 
implies the assumption of zero correlation of systematic 
errors among experiments (a reasonable first approxima- 
tion). 

To go beyond the first approximation, it is necessary to 
survey the existing experimental data on hadronic pro- 
duction cross sections. The published data and recent 
preprints on exclusive hadron production in electron- 
positron annihilations used in the present calculation are 
listed in Appendix ^ where it is evident that eight detec- 
tors have measured more than one mode. Since a given 
detector uses the same luminosity and similar correction 
factor calculations (e.g. radiative corrections, efficiencies 
which use some of the same subroutines) it is reasonable 
to suppose that the cross section determinations of differ- 
ent exclusive hadronic modes by a single detector may in 
fact be correlated. The following correlation postulates 
are therefore intended to address this situation: 

1. A single detector measuring more than one exclu- 
sive hadronic mode has 100 % correlations among 
systematic errors due to common luminosity and 
correction factor calculations; 



2. Different detectors have uncorrelated errors since 
they do not share luminosity and correction factor 
calculations. 

The accomodation of these correlation postulates re- 
quires a particular sequence of the three sums (1,2,3) for 
combination of the a^'"' central values and an alterna- 
tive sequence (1,3,2) for the a^'"' errors. In both cases, 
the energy is integrated over first and separately for each 
detector and exclusive hadronic mode measured over en- 
ergy sub-regions, where these sub-regions are defined by 
common energy coverage among detectors. 

The former sequence (1,2,3) facilitates the need to first 
calculate the error weighted fractions (defined in detail 
below) while the latter sequence (1,3,2) is necessary for, 
according to postulate 1, correlating individual detector 
systematic errors over the modes measured by that de- 
tector. This must be done before the final uncorrelated 
combination of errors is made, according to postulate 2, 
across detectors. As this method is based on Weighted 
Fraction averaging with S-factor application (see Eq. |^) 



After 
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integration over the energy, it is here referred 



to as the WFSA method. 

Lastly, in the present calculation it is noted that ex- 
clusive hadronic cross sections (modes) up to 2.0 GeV 
have been used because this can reveal the propagation 
of errors from each detector and exclusive mode sepa- 
rately. An additional consideration, although less impor- 
tant, is that the exclusive hadronic mode spectra may 
contain interference effects (when more than one vector 
meson contributes to a given exclusive hadronic mode) 
which would otherwise require special care if individual 
vector meson contributions were calculated separately. 
Further, the uncertainty over the generation mechanism, 
between the e"'"e~ annihilation photon and the hadronic 
final state, is avoided by focusing on exclusive hadronic 
modes themselves. 

For the energy region 2.0 to 3.1 GeV, the inclusive 
cross section ratio R{s) has been used (in the absence 
of exclusive data in this energy region) with the WFSA 
procedure for the contributing experiments. In the re- 
gion above 3.1 GeV, the QCD expression has been used 
without the WFSA procedure since p erturbative QCD is 
expected to be valid (see Sect. tVD). 



B. Trapezoidal Integration Procedure 

The usual trapezoidal integration technique takes the 
experimental data points pairwise: the cross sections, 
systematic and squared statistical errors are averaged per 
pair, then multiplied by the energy width of the pair, and 
finally they are summed over all pairs. However, it is con- 
venient for treatment of statistical errors to expand the 
sum in order to remove terms which cancel. Denoting by 
Sk, Kk, Cijk, cffjf, and cr^J^^ the energy, kernel function. 
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the cross section and its statistical and systematic error 
from the ith detector, jth exclusive mode, at the kth en- 
ergy point, the integration of Eq. (P can be represented 
by: 



1 V^'fl. ^U 

fc=l ^ ^ 



Sfe) 



k=2 



Ai = Cij^iKi{s2 - si) 

Ak = Cij^kKk{sk+i - Sfe-i) 

An — Cij^n-^n{Sn ^n—l^ 

where the first and last terms in the sum are handled sep- 
arately and the middle terms have an energy width across 
both upper and lower neighboring data points instead of 
across the points in pairs. This latter form is necessary 
for the proper treatment of the statistical errors: 
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ak'^iatl^^Kkfisk+i-Sk-.Y 

The systematic errors are specified by an array of val- 
ues p^j^l which are given as a percentage of the total cross 



section. Since the a,, contribution terms above are lin- 
ear in the cross section this implies that the systematic 
errors should be: 



sys sys a , \ ^ sys a , sys a 

% = P^llAl + 2^ P^lkAk + P^lnAn- 

k=2 

Usually the systematic error is the same for all energy 
points; however, it is different for each point of the NA7 
and Olya tt+tt^ measurements. In the cases where no 
systematic errors were given or readily located in the lit- 
erature, values from comparable measurements with the 
same detector have been taken wherever possible, or as 
10 % if the statistical errors dominate. In addition, a cal- 
culation with 20 % for these ambiguous values has shown 
that the WFSA results do not depend on the arbitrary 
selection of 10 % ||l^. All of the systematic errors that 
were used are listed in Tab. |l^ for tt+tt^, Tab. ^ for 



Tab. XI for the higher multiplicity modes, and 



Tab. XII for the energy region 2.0 - 3.1 GeV. 

If the desired limits of integration are inside (outside) 
the given data energy range, the cross section and error 
are linearly interpolated (extrapolated) to give the rele- 
vant pair of points. Thus our energy ranges are variable 
and have been set to match each of the previous evalua- 
tions for detailed comparisons p6[. 



C. Weighted Fraction Averaging 



In order to arrive at an a^"^'^ contribution per exclusive 

mode, a™"**^ = Oj, error- weighted averages across detec- 
tors have been performed according to the PDG ||l^ pro- 
cedure (for the ith detector and jth exclusive hadronic 
mode): 
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1/2 



where 



2 t ^stat\'^ , f^^y''\'^ 



(3) 



(4) 



In addition, the quality of the error weighted combina- 
tions were assessed by calculating the and PDG scale 
factor (i.e. per degree of freedom) |17t: 



N 
E 



N - 1 



(5) 



where N is the number of detectors included in the av- 
erage. If 5 > 1 then the errors were scaled up by this 
factor. 

At this point the prescription for determining the 
error-weighted fractional contribution to the total error 
from a given statistical or systematic error of the ith de- 
tector and jth exclusive mode is needed. For this purpose 
the PDG expression above for total squared error on a 
weighted average (see Eq. |[) can be expanded as a sum 
over squared component errors where in the first step the 
trivial sum Ei^i = is used to multiply by unity: 




The last step makes use of the definition Wij — l/f^j 
while the definition for the remaining afj in the numera- 
tor (see Eq. ||) leads to the separation of squared statis- 
tical from systematic terms in the sum: 
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Note that the PDG scale factor S has been inserted to 
emphasize the fact that the errors are to be increased if 
and only if the scale factor S > 1. (These expressions 
differ from those in the previous a^°'^ calculation by 
the factors 1/N and S.) 

Armed with these expressions it is possible to im- 
plement the correlation postulates for the ith detector 
by summing the averaged systematic errors (i.e. the 
weighted fractional systematic error contributions to the 
total errors in Eq. |^) linearly over contributing exclusive 
modes, while leaving weighted fractional statistical er- 
rors (Eq. ^ uncorrelated. The correlation postulates are 
simply executed by the following sums: 



the BCF collaboration (not apparently included pre- 
viously) has been divided by cr^^ to make them also R{s) 
values. 



B. Correlations Over Modes 

If the hadronic cross sections are measured over similar 
energy regions, as in the case of the > 2 hadrons multi- 
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yStat 



(8) 
(9) 



Note that itself is not directly used to determine the 
total WFSA error; to sum aj over j in quadrature would 
yield a total error which ignores correlations. This is 
what most previous a'^'^'^ calculations have done. 



IV. APPLICATION OF THE WFSA METHOD 

A. Correlations Over Energies 

In general, if different detectors measure hadronic cross 
sections over different energy regions, there must be dif- 
ferent error-weighted fractions and ^-factors for each 
common energy region (ith detector, jth energy region). 
Such a WFSA application over energies has been applied 
to the dominant a'*'"^ contributions from the hadronic 
modes: tt+tt^, tt+tt^tt'^, K'^K^, and Ki^Kg. An ex- 
ample of WFSA correlations over energy is shown in 



plicity cross sections shown in Tab. IV, then the WFSA 
method can simply be applied over the modes (ith de- 
tector, jth hadronic mode). In this case, a single error- 
weighted fraction and ^-factor suffices for evaluation of 
a^"''^ central value and errors. As before, the af}"''^ central 
values are error- weighted averaged across detectors (hor- 
izontally in the table) and summed over modes (verti- 
cally), while the detector specific statistical (systematic) 
errors are combined in quadrature (linearly) across, in 
this case, the modes (vertically). The errors are not to 
be combined across detectors until all modes have been 
treated separately and then taken together. It is this 
correlation over all modes, in particular correlating the 
TT+TT^ detector total systematic errors with all the other 
two body and higher multiplicity modes, over all energies, 
which appears not to have been done before. (The pre- 
vious a^'^'^ calculation ||l^ appears to have accounted for 
these correlations only above 0.81 GeV, after the peak 
of the p-meson. Hence, it appears that the tt+tt" sys- 
tematic errors below 0.81 GeV were not correlated with 
systematic errors of measurements by the same detectors 
of the other modes over all energies and of the tt+tt^ 
mode above 0.81 GeV.) 



C. Kaons, Narrow Resonances and Higher VMD 
modes 

Some additional hadronic modes which contribute to 



Tab. ^ for the R{s) a'^'"^ contribution in the energy range the total hadronic production cross section CT°(e 



from 2.0 to 3.1 GeV. The a'},"'^ central values are error- 
weighted averaged across detectors (horizontally in the 
table) and summed over energies (vertically) while the 
detector specific statistical (systematic) errors are com- 
bined in quadrature (linearly) across the energies (verti- 
cally). The correlation of systematic errors over energies 
does not appear to be taken into account by the WFB 
method discussed in the previous ajj"'* calculation p2[ : 
Weighted fraction averaging (of R{s) values) Before en- 
ergy integration. 

In this energy range, published inclusive R{s) val- 
ues have been used from the 772 detector Q jl^] and 
MARK I [|0| . It is noted that the 772 collaboration has 
published values for i?2(s) (two hadrons exclusively) |l^] 
and i?>3(s) (three or more hadrons) jl^] separately which 
are added for the present calculation. This does not ap- 
pear to have been done in the previous a'^"'' calculations. 



In addition, data published as (T(e^ 



hadrons) by 



hadrons) require comment. In particular, experimental 
data in the form of total cross sections neither exists on 
the radiative decays of the ui and (j) mesons (7r'^7, 777) nor 
on the kaon pair production of the (p meson below certain 
energies. Further, there are additional contributions to 
hadronic vacuum polarization than from just the lowest 
order single vector meson dominated (VMD) amplitudes 
represented by the hadronic (decay) modes previously 
discussed. 



1. Kaons and Narrow Resonances 

Experimental data on the total cross sections for pro- 
duction of kaon pairs (charged and neutral) are limited 
by the fact that nuclear interactions of low momentum 
kaons are not well measured. Hence kaon detection effi- 
ciencies are difficult to calculate precisely, and although 
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existing data does span both sides of the </) meson it is 
usually a relative cross section useful for measurements 
of (/) meson parameters. Hence a Breit-Wigner (BW) line 
shape with PDG 1994 parameters has been used, below 
the (most recent) lowest data point for kaon pair pro- 
duction and for the whole energy range for the radiative 
decays, where the errors (considered totally systematic) 
are evaluated by differentiating the BW formula. 

In the charm and bottom threshold regions the 6 states 
each of the J/^ and T resonance families have been cal- 
culated separately using PDG 1994 values and the peak 
approximation formula for a'^"''^ ( ||2|] ||ll|). In view of 
the small contribution to a|j'"^ from the &-quark (due to 
kernel function suprcssion in Eq. |l]), the top quark con- 
tribution is neglected. 



2. Higher Order VMD Contributions 



The Vector Meson Dominance (VMD) model approx- 
imation of QCD expresses the fact that vector mesons, 
instead of quarks and gluons, are the relevant degrees 
of freedom in low energy (near threshold) QCD inter- 
actions. The VMD model approximation for hadronic 
vacuum polarization is depicted in Fig. |[ Most of the 
^had contributions arc contained in the first term on the 
right of the figure. The decay modes of the vector mesons 
p, CO, <j) and ui' account for - 597 x lO^^" (85% of a'/,'"^), 
while p' accounts for - 34.7 x 10^^° (5% of a'/^"'^), and 
(j)' for ^ 1.6 X 10^"'^° by rough accounting. The remain- 
ing part (~ 10%) of ajj°'^ is derived from a non- or very 
broadly resonant background for which it is instructive to 
consider some additional generation mechanisms not al- 
ready included in the previous exclusive mode approach. 

Namely, some portion of a^"'* derives from higher or- 
der VMD interactions Vtt and Vtttt (the next terms on 
the right in Fig. ||) where F is a vector meson. Although 
these amounts are individually less than the expected 
E821 error (~ 4.0 x lO"^*^) the fact that they may com- 
bine with other small contributions (whose sum may be 
greater than the E821 error) implies that they all should 
indeed be carefully considered. 

The cross sections a{K* K^tt^), cr(w7r''), a{ujTr~^Tr~), 
cr(?77r+7r~) have recently (1991,1992) been measured 
(ND,DM2). In particular, it has been pointed out p^ ] 
that the following exclusive hadronic modes marked with 
the comment " not " need to be included (B refers to the 
branching fraction for the decay mode in parentheses): 
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is in TT+TT" TT^TT^. 



They further point out that isospin considerations |23| ] 
imply 

2ct(tT+TT-Tt"tt"tT°) = Cr(TT+TT-TT+TT-TT") 
2cr(a'TT"TT°) = fT(wTT+TT-) 

and that this must be used since ct(wtt*'tt*') and 
(t(tt+tt-tt°tt°tt'') have not been measured. Therefore the 
cr(TT+TT-TT+TT-TT'^) Contribution has been augmented by a 
factor Ai = 1.5, and the contributions from a{ujT:^) and 
cr(a;TT+TT-) will contribute with factors Bi = 0.085 and 
= 1.5Bi respectively, as noted in Tab. [V. However, 



the same logic applied to cr(?7TT+TT-) implies not including 
100% of it as done previously [l^ : 

(77)tt+tt- B=0.388 not in tt+tt^ 



TT+TT 



(tt"tt"tt")tt+tt- B=0.319 is in tt+tt-tt"ttOtt° 
(TT+TT-TT'^)TT+Tr- B=0.236 is in tt+tt^tt+tt^tt*^ 



(tt+tt 7)tt+tt B=0.0488 not in tt+tt tt+tt tt°. 
To avoid double counting in the tt+tt-tt+tt^tt" and 
tt+tt-tt'^tt^tt'^ channels only the fraction B2 — 0.388 -I- 
0.0488 of the o'(?7tt+tt-) contribution has been included. 
(There is no further factor of I /2 for inclusion of the cross 
section a{e~^e~ — > TyTT^Tr") since it is forbidden.) The 
augmentation and branching factors {Ai, Bi, B2, B3) 
are noted in Tab. ^ where all contributions are listed in 
descending order. 



D. The Perturbative QCD Energy Region 



To test the validity of QCD for determination of the 
contribution to a^/f-'^ in the energy region above 3.1 GeV, 
a QCD parameterization [|4j (including second order 
terms) was compared with a data-based evaluation p2[ . 
An asymptotic kernel function ||l^ was used in the a^'^ 
integral since Eq. ^ is numerically stable only up to ~ 20 
GeV. The errors on the QCD contribution to a'^ " ' 



^had 



determined by its variation with A-p^ ± AA^^^^. 

Marshall combined the results from 15 different e+e" 
annihilation experiments, fitting them to a third order 
QCD model with a single parameter, A-^g-j. Because the 
fit was overconstrained, he was able to evaluate whether 
these experiments had overestimated their systematic er- 
rors. He then went a step further, and fitted for the 
absolute normalization of each of the 15 experiments in- 
dependently, bounded by double their stated systematic 
errors (±2(t). He found fitted normalizations for most ex- 
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periments within their stated hmits, with the exception 
of two of the earhest experiments: Mark II and 772. 

By using the error in his normahzation constants (from 
the fit) rather than the stated systematic errors for each 
experiment, Marshall was able to significantly reduce his 
error on Ajjg- and hence on his overall normalization for 
R{s). Dubnicka followed Marshall in his 1992 preprint 
1^, and in the present calculation both Marshall's fitted 
parameterization (with his errors) and trapezoidal inte- 
gration have been compared for the higher-energy 
contributions to ajj"''. As the QCD and data-based re- 
sults are in good agreement it is clear that the second 
order QCD expression is sufficient. 



V. WFSA RESULTS AND PROJECTIONS 



the neglect of third order QCD terms; as new data bring 
the errors from the lower energy regions downward (see 
Sect. 0), the QCD error will become more significant 
however.) 

It is further apparent from Table ^ that the error in the 
energy region below 1.4 GeV dominates the total error in 



-.had 



and this error comes mostly from the Olya tt^ 



data. (One advantage of the WFSA method is that in 
presenting averaged errors (both statistical and system- 
atic separately) for each detector and exclusive mode, it 
remains possible to identify which detector and exclusive 
mode is contributing the most to the overall total error.) 
The WFSA method most heavily weights the data with 
the smallest errors (see Table IX ) , and hence most of the 
stated systematic errors are indeed being added linearly 
for the experiments with the smallest errors. 



A. Summary Results 

As the preceding discussion of the WFSA evaluation 
shows, the core ajj"*^ problem is how to combine the cor- 
relations among the systematic errors. In the present cal- 
culation, partial a^"'^ integrations have been performed 
first so relative errors and error-weighted averaged sys- 
tematic and statistical errors could be defined at an early 
stage of the error combination calculation. The aver- 
aged systematic errors have bee n corr elated according to 
the postulates discussed in Sect. Ill A for all experiments 
measuring more than one hadronic mode; all other er- 
rors have been combined in quadrature only after such 
correlations have been taken into account. 

All of the WFSA calculations have been collected in 
Tab. ||. The first two lines of the table are summary 
results of the a^"'' contributions and errors presented in 



Table |l| In the > 2 hadrons line of Tab. |l| the results 
from Tab. IV are presented. The third line line in Tab. || 
for R(s) from 2.0 - 3.1 GeV is the summary result from 

Table All of the aj^"'* central values are averaged over exclusive hadronic modes: 



B. Energy Region < 1.4 GeV 

The overall WFSA results indicate that improved mea- 
surements below 1.4 GeV, can significantly improve upon 
the overall lowest order hadronic vacuum polarization 
uncertainty in a'^"''^. To show this, the WFSA proce- 
dure has been performed with an additional two ex- 
periments in view, the CMD2 experiment at VEPP-2M 
in Novosibirsk, Russia and a hypothetical detector at 
DAFNE at Frascati, Italy. The a'^'"^ resuhs for the en- 
ergy region below 1.4 GeV are shown in Tab. VI for 
three cases: a) without new data, b) with expected er- 
rors (stat)(sys) = (3%)(0.5%) from the CMD2 experi- 
ment and c) CMD2 plus a second experiment with errors 
(stat)(sys) = (3%)(0.5%). In particular, the central val- 
ues for hadronic cross sections have been chosen equal 
to CMD or Olya values (and hence the S factors will 
generally be less than one and so not applied) , while the 
errors above have been determined for only the following 



(2 tt) 



(3 tt), 



detectors first and then summed over modes, while the 
errors are combined over modes first (to correlate system- 
atics) and only combined over detectors in the end. The 
errors m Tab. may be combined in quadrature since 
there are no correlations remaining among the categories 
chosen. (This feature is not completely present in the 
previous calculations of a^""^.) 

A subtotal for aj^'"' contributions below 3.1 GeV is pre- 
sented in Tab. || in order to show the total errors before 
the use of QCD. It is there seen that the QCD results 
do not have much influence on the final results since the 
error from the lower energy (< 3.1 GeV) region domi- 
nates. For comparison, if in place of the QCD calcula- 
tion the data based evaluation for the region above 3.1 
GeV [|l2| is used the final errors are only slightly higher: 
^had ^ 702.718(7.787)(14.147). Hence, the WFSA result 
is thus shown to be largely insensitive to the use of QCD 
or the experimental data after 3.1 GeV. (This justifies 
both the use of QCD down to the energy 3.1 GeV and 



(4 tt), TT+TT TT TT (4 7r), TT+TT TT+TT TT (Stt), 

K+K- and KlKs. 

The choice of 3 % statistical error is based on the 
CMD2 data taking assumption of 1000 pion pairs per 
energy point since tt+tt^ is the dominant contribution 
to ajj°'^. Since the luminosity collected by CMD2 has 
been determined by this requirement on the tt~^7t~ cross 
section (and by competing requirements for other vec- 
tor meson physics goals), the actual statistical errors on 
the non-TT+TT^ higher multiplicity modes will in fact be 
somewhat different, but this is a higher order effect here 
neglected. 

The choice of 0.5 % systematic error is based on the 
fact that the limiting error on the new CMD2 cross sec- 
tion measurements appears to be the error on higher or- 
der corrections to the QED Bhabha cross sections used 
in calculations of the luminosity. While radiative correc- 
tions have been calculated to good accuracy (0.11 % [p5[) 
for the i-channel contributions to Bhabha scattering use- 
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ful for the forward region luminosity monitors in use at 
LEP, they are not useful for CMD2; the principle lumi- 
nosity is determined there by large angle Bhabha events 
in the barrel calorimeter where significant t- and s- chan- 
nel interference terms are present. The task of calculating 
these interference terms to 0.5 % accuracy is well under- 
way and it is assumed that this will be the limiting 
error on the new CMD2 cross section measurements. 

The results in Tab. ^ show that without new data the 
errors on the contributions to 10^°a'},'"^ below 1.4 GeV 
will be 10^° (stat)(sys) = (7.4)(13.1) and with particular 
new data they can become 10^° (stat)(sys) = (2.1)(2.9) 
which is equivalent to a total error of 3.6 x lO"""^". Com- 
paring this with the AGS E821 experimental error goal 
of (Texp ^ ± 4.0 X 10~^° it is clear that the new data is 
needed. In particular, if this data has errors of 3 % sta- 
tistical and 0.5 % systematic on the stated modes, then 
they will in fact be sufficient to reduce the hadronic con- 
tribution errors in the energy range below 1.4 GeV to 
below the error goal of the new measurement of ajj°'^ by 
the BNL AGS E821 Collaboration. 



C. Energy Region > 1.4 GeV 



New data on the 4 tt neutral, 5 tt and 6 tt modes can 
reduce the 772 contribution to the total errors (as shown 
in Row 2) while new data on the 4 tt charged mode is 
required to reduce the DM1 and DM2 contribution to 
the total errors (as shown in Row 3). Therefore, if many 
data points 10 — 20) for each of the specified modes 
across the energy region 1.4 to 2.0 GeV with the assumed 
errors can be extracted from r-decay data, then the con- 
tribution to the total errors from this ene rgy region will 
be reduced significantly as shown in Tab. VIl. 



2. New R{s) Measurements for Energy Region 2.0 - 3.1 GeV 



The effect of new data on the error in the 2.0 - 3.1 GeV 



energy region is presented in Tab. VIIl. The WFSA pro- 
cedure has been performed with a new hypothetical de- 
tector measuring R{s) from 2.2 - 3.1 GeV (central value 
of R[s) — 3.0 for all points so no ^-factors exceed 1) with 
assumed errors of (stat)(sys) = (3%)(0.5%). If new R[s) 
data in this energy region can be obtained with these er- 
rors then the 772 contribution to the total error can be 
significantly reduced. 



However, the new data from BINP Novosibirsk and 
INFN Frascati will not reduce the errors on a'^"''^ in the 

energy range above 1.4 GeV. The current error on ajj'"' 
obtained by use of the WFSA procedure in the region 
above 1.4 GeV is lO^" (stat)(sys) = (2.4)(5.3). In this 
region, most of the error comes from the 772 detector 
at the former 3 GeV Adone storage ring in INFN Fras- 
cati, Italy. As this error is larger than the BNL AGS 
E821 experimental error goal of atxp ± 4.0 x 10"^", 
clearly more experiments are needed for interpretation of 
the new measurement of ajj:'"^. Fortunately, it is possible 
to use new high statistics data on r-decays for new mea- 
surements of multi-pion production in the energy region 
1.4 - 2.0 GeV, and above 2.2 GeV there are plans to make 
further measurements of the cross section ratio R{s) at 
BEPC Beijing, China 0. 



1. T Decay Data for Energy Region I.4 - 2.0 GeV 

The new high statistics data on r-decays has already 
been included in the present calculation by use of the 
bjn^ cross section obtained from the ARGUS Collabo- 
ration 1^ ll^]. The idea is based on the fact that the 
coupling of bosons to quarks is related to the photon 
coupling to quarks by an isospin rotation (CVC relation) 
pO| . The effect on the total errors in the energy region 
1.4 - 2.0 GeV are presented in Tab. VII where the as- 
sumed errors are (stat)(sys) = (3%)(1%) for the following 
hadronic modes: 7r+7r^7r+7r^ (4 tt charged), 1:^1:^1:^1:^ 



D. Conclusion 



Taking the new hypothetical data together (including r 
decays and R{s) measurements), the total error from the 
energy range above 1.4 GeV (including resonances and 
QCD) then becomes 10^" (stat)(sys) = (0.6)(1.6). Com- 
bining the hypothetical results from above and below 1.4 
GeV the total error becomes 10^° (stat) (sys) = (2.2) (3.3) 
which is equivalent to a total error of 3.98 x 10~^°. There- 
fore, new measurements of e^e~ hadrons (including 
the cross section ratio R{s)) and r decays with assumed 
errors of (stat)(sys) = (3%)(0.5%) and (3%)(1.0%) re- 
spectively are sufficient to reduce the total error on the 
lowest order contribution to aj^"'' in all energy regions 
below the expected error of the new measurement of 
by the AGS E821 CoUaboration. 



(4 TT neutral), tt'^tt tt 
(6 tt). 



(5 tt) and 1:^1: tt+tt tt tt' 
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APPENDIX A: EXPERIMENTAL DATA USED 



The experimental data we used for our evalua tion of 

-,had • 



is listed by hadronic exclusive mode in Tab. XIII 
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TABLE I. History of aj;"'* calculations. The values 4 and 15 x 10"^° are the AGS E821 error goal and electroweak contribution 
to respectively. 



Calculations based on Experimental Data 


Author(s) 


Year 




ref 


Budker Institute 
Kinoshita, Nizic, Okamoto 
Dubnickova, Dubnicka, Stricinec 
Eidelman, Jegerlehner 
Apel, Yndurain 
WFSA evaluation 


1985 
1985 
1992 
1995 
1995 
1996 


684(11) 
707(6) (17) 
699(4X2)" 
702(6) (14) 

710(11) 
703(8) (14) 







Improvement due mostly to fitting function technique. See text. 



TABLE IL The Grand Total WFSA lO^^af^'"' resuhs 



Energy Region [GeV] 


WFSA lO^^aJi"" (stat) (sys) 


(7(e"''e~ hadrona) Ee - 1.4 
(j(e^e^ hadrons) 1.4 - 2.0 

R(s) 2.0 - 3.1 

J/* (6 states) 
T (6 states) 

QCD 3.1 - oo 


611.332 (7.399) (13.045) 
32.466 (0.756) ( 2.379) 
28.374 (2.288) (4.400)" 
9.047 (-) (0.969) 
0.109 (-) (0.013) 
21.301 (-) (0.371)'' 


Sub Total < 3.1 + J/*,T 


681.328 (7.782) (14.005) 


Sub Total < 1.4 


611.332 (7.399) (13.045) 


Sub Total > 1.4 


91.297 (2.410) ( 5.108) 


Total Figure 1 


702.629 (7.782) (14.009) 



Eg refers to tlie particular thresholds of the exclusive hadronic modes. 

° Represents systematic errors mainly from 772 detector added linearly for energy regions 1.4 - 2.0 and 2.0 - 3.1 GeV 
(1.7 + 2.7). 

^ Errors determined by a'^'^^ {^TTs ^ ^^TTs) 
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TABLE III. Our WFSA lO^^a^'"' results below 2.0 GeV. Upper (lower) contributions listed in the total column are for below 
(above) 1.4 GeV. Upper (lower) numbers in parentheses are contributions to the statistical (systematic) errors. Systematic 
(statistical) errors are combined linearly (in quadrature) in each column separately. 



Mode 


CMD 


Olya 


ND 


DMIA 


DMID 


DM2 


Other 


Total 




(4.806) 
(4.789) 


(1.306) 
(5.807) 




(5.064) 
(4.021) 




(0.058) 
(0.110) 


(0.814) 
(0.732) 


502.184 
0.781 




(0.509) 
(0.104) 




(0.856) 
(1.303) 


(1.209) 
(0.715) 




(0.041) 
(0.053) 


(0.049) 
(0.009) 


50.886 
0.677 


> 2 hadrons 


(0.236) 
(0.589) 


(0.369) 
(2.405) 


(0.205) 
(2.070) 




(0.241) 
(1.071) 


(0.307) 
(1.857) 


(0.631) 
(1.735) 


16.782 
29.995 


K+R- 


(0.119) 
(0.058) 


(0.953) 
(1.900) 








(0.046) 
(0.090) 


(0.230)'' 


20.623 
4.45" 
0.759 


KlKs 


(0.173) 
(0.040) 


(0.156) 
(0.101) 






(0.033) 
(0.017) 




(0.586)'' 


0.755 
14.07" 
0.154 


PP 










(0.020) 
(0.010) 






0.100 


ui 7r"7, rij 














(0.040)" 


0.980" 


cj) —> 7r^7, 777 














(0.024)" 


0.602" 


Total 
Ee - 1.4 


(4.842) 
(5.580) 


(1.666) 
(10.213) 


(0.880) 
(3.373) 


(5.206) 
(4.736) 






(0.807) 
(0.962) 


611.33 
(7.399) 
(13.045) 


Total 
1.4 - 2.0 










(0.244) 
(1.098) 


(0.318) 
(2.110) 


(0.642) 
(1.747)^= 


32.47 
(0.756) 
(2.379) 



Eg refers to the particular thresholds of the exclusive hadronic modes. 
DMIA = DM1 at AGO; DMID = DM1 at DCI, Orsay. 



" Integration of energy dependent width Breit-Wigner in absence of data. 
^ Errors determined by my, Fy, By^ee derivatives of the Breit-Wigner. 

^ The 772 systematic errors from energy regions 1.4 - 2.0 GeV (1.7) and 2.0 - 3.1 GeV (2.7) will be added linearly 
and presented in Table |l[ 
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TABLE IV. Our WFSA lO^^aJj"'' results for > 2 hadrons exclusive modes below 2.0 GeV. Upper (lower) contributions listed 
in the total column are for below (above) 1.4 GeV. Upper (lower) numbers in parentheses are contributions to the statistical 
(systematic) errors. Systematic (statistical) errors are combined linearly (in quadrature) in each column separately. 



M()(l(- 


CMD 


Olya 


Nl) 


DA Hi) 


1)A12 


Otlit^r 


Total 


S Factor 


7r"'"7r~7r"7r" 




(0.325) 
(1.491) 


(0.193) 
(1.514) 




(0.256) 
(0.868) 


(0.335) 
(0.841) 


10.783 
8.543 


1.614 
1.692 




(0.225) 
(0.527) 


(0.174) 
(0.914) 


(0.066) 
(0.526) 


(0.104) 
(0.623) 


(0.070) 
(0.628) 




5.161 
10.215 


0.486 
1.221 


7r"'"7r~7r"'"7r~7r"7r" 












(0.494) 
(0.763) 


5.089 


0.0 


^i(7r+7r~7r+7r^7r'-') 


(0.072) 
(0.062) 






(0.124) 
(0.192) 




(0.187) 
(0.131) 


0.305 
2.533 


0.0 
1.010 










(0.158) 
(0.0) 


(0.104) 
(0.119) 




0.951 


3.199 


K+K-TT+iv- 








(0.072) 
(0.135) 


(0.091) 
(0.123) 




0.815 


2.709 












(0.057) 
(0.069) 




0.692 


0.0 








(0.017) 
(0.030) 






(0.079) 
(0.021) 


0.533 
0.210 


1.132 
0.0 












(0.039) 
(0.044) 




0.444 


0.0 


7r"'"7r~7r"'"7r~7r"'"7r^ ) 








(0.040) 
(0.117) 






0.419 


0.0 










(0.005) 
(0.003) 


(0.003) 
(0.004) 




0.084 


0.982 


Total 


(0.236) 
(0.589) 


(0.369) 
(2.405) 


(0.205) 
(2.070) 


(0.241) 
(1.071) 


(0.307) 
(1.857) 


(0.631) 
(1.735) 


16.782 
29.995 





TABLE V. Our WFSA 10^°aJ;'"' results for R{s) from 2.0 - 3.1 GeV. Error weighted averages with S-factor in "Total" column. 

Numbers in upper (lower) parenthesis are WFSA contributions to the statistical (systematic) errors which are combined in 
quadrature (linearly) in each column separately. Lower right corner errors are quadrature sums of errors in the "Total" row. 



Energy Range [GV] 


EOF 


772 


Mark I 


Total 


S Factor 


2.0 - 2.6 


18.778 
(2.149) 
(0.207) 


22.441 
(0.527) 
(2.093) 




20.480 


0.637 


2.6 - 2.87 


3.718 

(0.469) 
(().():!2) 


6.535 

(0.130) 
(0.151) 


4.781 

(0.2.36) 
(O.IDG) 


4.898 


1.210 


2.87 - 3.0 


1.721 
(0.202) 
(0.019) 




1.777 
(0.061) 
(0.193) 


1.742 


0.096 


3.0 - 3.1 






1.254 
(0.062) 
(0.313) 


1.253 


0.0 


Total 


(2.208) 
(0.258) 


(0.543) 
(2.544) 


(0.251) 
(0.913) 


28.374 
(2.288) 
(2.700) 
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TABLE VI. WFSA lO^^aJj"'' detector total statistical (upper) and systematic (lower) errors below 1.4 GeV: Row 1) with no 
new data, Row 2) one new experiment with 3 % stat., 0.5 % syst. errors, Row 3) two new experiments with 3 % stat. and 0.5 
% syst. errors. 









WFSA 10'%;;;" 


" Errors Below 1.4 GeV 






iNew z 


New i 




Ulya 


iNJJ 


UMiA 


Utner 


iotal 






(4.842) 


(1.666) 


(0.880) 


(5.206) 


(0.807) 


(7.399) 






(5.580) 


(10.213) 


(3.373) 


(4.736) 


(0.962) 


(13.045) 




(2.426) 


(1.266) 


(1.148) 


(0.136) 


(1.239) 


(0.451) 


(3.250) 




(2.224) 


(1.292) 


(3.608) 


(0.405) 


(1.241) 


(0.760) 


(4.681) 


(1.183) 


(1.183) 


(0.789) 


(0.528) 


(0.078) 


(0.766) 


(0.361) 


(2.103) 


(1.256) 


(1.256) 


(0.773) 


(1.921) 


(0.244) 


(0.769) 


(0.718) 


(2.934) 



TABLE VII. WFSA 10^°oji'"' detector total statistical (upper) and systematic (lower) errors in region 1.4 - 2.0 GeV. Row 1) 
with no new data. Row 2) new hypothetical data (possibly from t decay spectral functions) on modes 7r"'"7r~7r°7r°, 7r"'"7r~7r"'"7r~7r° 
and 7r"'"7r~7r"'"7r~7r''7r° with 3 % stat., 1.0 % syst. errors. Row 3) including as well the mode 7r"'"7r~7r"'"7r~. 







WFSA 10'"a;:'"' Errors 


in 1.4 - 2.0 GeV Region 




New 1 


DMID 


DM2 


772 


Total 




(0.244) 


(0.318) 


(0.628) 


(0.745) 




(1.098) 


(2.110) 


(1.744) 


(2.950) 


(0.068) 


(0.210) 


(0.192) 


(0.067) 


(0.300) 


(0.115) 


(0.922) 


(1.320) 


(0.141) 


(1.620) 


(0.082) 


(0.179) 


(0.157) 


(0.067) 


(0.261) 


(0.175) 


(0.346) 


(0.514) 


(0.141) 


(0.659) 



TABLE VIII. WFSA lO^aJ^'"^ detector total statistical (upper) and systematic (lower) errors in region 2.0 - 3.1 GeV. Row 
1) with no new data. Row 2) one new experiment above 2.2 GeV with 3 % stat., 0.5 % syst. errors 







WFSA lO'^o^l'"' Errors 


in 2.0 - 3.1 GeV Region 




New 1 


BCF 


772 


Mark I 


Total 




(2.208) 


(0.543) 


(0.251) 


(2.288) 




(0.258) 


(2.544) 


(0.913) 


(2.715) 


(0.109) 


(0.485) 


(0.140) 


(0.025) 


(0.517) 


(0.058) 


(0.109) 


(0.601) 


(0.084) 


(0.619) 



14 



TABLE IX. S^'Htcmatic Errors used in our evaluation of the tt tt coutriliutioii to a ' 



Detector 



Systematic error 



NA7 
CMD 
TOF 
DM1 

Olya 
M2N 
BCF 

fill 
MEA 
D_\I2 



.01,.02,.025,.05 
.02 
.035" 
.022 
.04 - .15 
.0" 
.1= 
.1= 

.r 

.12 



" evaluated from errors on radiative eorreetion/effieiency factors. 

^ systematic error is included with statistical in quoted errors. 

" Arbitrary value - inconsequential since statistical errors are large. 

TABLE X. Systematic Errors used in our evaluation of the 'k~^'k~tP contribution to a^"'*. 



Detector 


Energy Range [GeV] 


Systematic error 


CMD 


.76 - .81 


.0" 


CMD 


.84 - 1.013 


.07 


ND 


.414 - .765 


.1 


ND 


.805 - 1.003 


.1 


ND 


1.036 - 1.379 


.2 


DM1 


0.414 - 1.098 


.032 


DM2 


1.34 - 2.0 


.0866 


772 


1.437 - 2.0 


.15 



" systematic errors are less than quoted statistical errors. 
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TABLE XI. Systematic errors used in our evaluation of higher multiplicity exclusive hadronic mode ajf contribution. 



Mode 


Olya 


CMD 


ND 


ARGUS 


DM1 


DM2 


772 




o.r 


0.038 


0.1 






0.1" 




KlKs 


0.3 


0.1 


0.1 




o.r 








0.15 


0.1 


0.1 




0.1 


0.1" 




+ -00 
TT TT IT TT 


0.15 




0.15 






0.1 


0.15 


+ - + - 
TT 71" n TV TT 




0.20 






0.12 




0.15 


K+K-TT+TT- 










0.13 


0.1" 














0.0 


0.1" 
















0.1" 














0.28 






+ - + -00 

n n TV n n n 














0.15 


PP 










0.1" 






u)n 






0.1 


0.1 


















0.12 


0.082 
















0.1" 





" systematic error not discussed in reference 



TABLE XII. Systematic errors used in evaluation of the energy region 2.0 - 3.1 GeV a^""* contribution. 



Energy Region [GeV] 


BCF 


772 


Mark I 


2.0 - 3.1 


0.02 


0.21 


0.25 
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TABLE XIII. Data used in the WFSA evaluation 



Mode 


OLYA 


CMD 


ND 


DM1 


DM2 


Others 


t: 

■K TV 


H 


i 






PI 


M2N BCF ^ivr MEA NA7 TOE 

H N H H H H 


TT TT TT 




|40l P 








GG281 














MEA 
PI 


KlKs 














T X 

TT 71" TT TT 




|49| 


PI 




1^ 




X n — n 


PI 




PI 






GG281 

m 














GG281 


T T T 

TT TT TT TT TT TT 








[R 






T TT n — n — 

TT TT TT TT TT TT 












GG281 












1^ 




mr^TT X 

K+K TT + TT 
























¥% 










PI 






ARGUS 
PI 


t: 










a 




CJTTTT 








PI 


PI 


















i?(s) 2.0 - 3.1 GeV 












GG279 GG281 BCF74 MarkI 

@ @ ill 
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FIG. 1. Lowest order hadronic contribution to aj^"''. 




